
# Read the set file
raw_path = './oddball/sub-001/eeg'
raw = mne.io.read_raw_eeglab(raw_path + "/sub-001_task-P300_run-1_eeg.set", preload = True)
raw


# Visualize the raw data
%matplotlib inline
raw.plot();


# Number of channels before subsetting
len(raw.ch_names)

79


# Show all channels before subsetting
', '.join(raw.ch_names)

'Fp1, AF7, AF3, F1, F3, F5, F7, FT7, FC5, FC3, FC1, C1, C3, C5, T7, TP7, CP5, CP3, CP1, P1, P3, P5, P7, P9, PO7, PO3, O1, Iz, Oz, POz, Pz, CPz, Fpz, Fp2, AF8, AF4, AFz, Fz, F2, F4, F6, F8, FT8, FC6, FC4, FC2, FCz, Cz, C2, C4, C6, T8, TP8, CP6, CP4, CP2, P2, P4, P6, P8, P10, PO8, PO4, O2, EXG1, EXG2, EXG3, EXG4, EXG5, EXG6, EXG7, EXG8, GSR1, GSR2, Erg1, Erg2, Resp, Plet, Temp'


# Select only a subset of the channels
raw_eeg = raw.copy()
raw_eeg.pick_types(meg = False, eeg = True, eog = True, exclude = [])

NOTE: pick_types() is a legacy function. New code should use inst.pick(...).


# Number of channels after subsetting
len(raw_eeg.ch_names)

65


# Show the subset of the channels
', '.join(raw_eeg.ch_names)

'Fp1, AF7, AF3, F1, F3, F5, F7, FT7, FC5, FC3, FC1, C1, C3, C5, T7, TP7, CP5, CP3, CP1, P1, P3, P5, P7, P9, PO7, PO3, O1, Iz, Oz, POz, Pz, CPz, Fpz, Fp2, AF8, AF4, AFz, Fz, F2, F4, F6, F8, FT8, FC6, FC4, FC2, FCz, Cz, C2, C4, C6, T8, TP8, CP6, CP4, CP2, P2, P4, P6, P8, P10, PO8, PO4, O2, Temp'


# In 2D
%matplotlib inline
mne.viz.plot_sensors(raw.info, kind = 'topomap', ch_type = 'all', show_names = True);


# Referencing using the T7, T8 channels as the referenced electrode
raw_ref = raw_eeg.copy()
raw_ref.set_eeg_reference(ref_channels = ['T7', 'T8'])

EEG channel type selected for re-referencing
Applying a custom ('EEG',) reference.


raw_ref.plot();


# Plot the power spectrum before filtering
raw_ref.plot_psd();

NOTE: plot_psd() is a legacy function. New code should use .compute_psd().plot().
Effective window size : 8.000 (s)


# Using a band pass filter from 0.16-30 Hz
raw_filtered = raw_ref.copy()
raw_filtered = raw_filtered.filter(l_freq = 0.16, h_freq = 30)

Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 0.16 - 30 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 0.16
- Lower transition bandwidth: 0.16 Hz (-6 dB cutoff frequency: 0.08 Hz)
- Upper passband edge: 30.00 Hz
- Upper transition bandwidth: 7.50 Hz (-6 dB cutoff frequency: 33.75 Hz)
- Filter length: 5281 samples (20.629 s)

[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done   1 out of   1 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   2 out of   2 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   3 out of   3 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   4 out of   4 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done  65 out of  65 | elapsed:    0.2s finished


fig, ax = plt.subplots(2, figsize = (10, 5))

# Plot PSD before filtering
raw_ref.plot_psd(ax = ax[0], picks = ['eeg'], show = False)
ax[0].set_title('PSD before filtering')

# Plot PSD after filtering
raw_filtered.plot_psd(ax = ax[1], picks = ['eeg'], show = False)
ax[1].set_title('PSD after filtering')

# Set x-axis label
ax[1].set_xlabel('Frequency (Hz)')

# Adjust layout and show the figure
fig.tight_layout()
plt.show()

NOTE: plot_psd() is a legacy function. New code should use .compute_psd().plot().
Effective window size : 8.000 (s)
NOTE: plot_psd() is a legacy function. New code should use .compute_psd().plot().
Effective window size : 8.000 (s)


%matplotlib inline
raw_filtered.plot_psd(fmax = 50);

NOTE: plot_psd() is a legacy function. New code should use .compute_psd().plot().
Effective window size : 8.000 (s)


# Get events and event ids from annotations
events, event_id = mne.events_from_annotations(raw_filtered)

Used Annotations descriptions: ['ignore', 'noise', 'noise_with_reponse', 'oddball', 'oddball_with_reponse', 'response', 'standard']


# Rejecting epochs with peak-to-peak amplitude +/-100 µV
reject_criteria = dict(
    eeg = 100e-6,  # 100 µV
)


# Epochs will start 200 ms before the stimulus onset and end 
# 600 ms after stimulus onset
tmin = -0.2
tmax = 0.6
baseline = (None, 0)

epochs = mne.Epochs(raw_filtered,
                    events = events,
                    event_id = event_id,
                    reject = reject_criteria,
                    tmin = tmin,
                    tmax = tmax,
                    baseline = baseline,
                    preload = True)
epochs

Not setting metadata
863 matching events found
Setting baseline interval to [-0.19921875, 0.0] s
Applying baseline correction (mode: mean)
0 projection items activated
Using data from preloaded Raw for 863 events and 206 original time points ...
    Rejecting  epoch based on EEG : ['P10']
    Rejecting  epoch based on EEG : ['P10']
    Rejecting  epoch based on EEG : ['FT7', 'P10']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FT7', 'FC5', 'FC3', 'FC1', 'C1', 'C3', 'C5', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC3', 'FC1', 'C1', 'C3', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC5', 'FC3', 'FC1', 'C1', 'C3', 'T7', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'T8', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC3', 'FC1', 'C1', 'C3', 'T7', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'T8', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC1', 'C1', 'T7', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'T8', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'F5', 'T7', 'TP7', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'C2', 'C4', 'C6', 'T8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'T7', 'P9', 'O1', 'Iz', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F8', 'FT8', 'FC6', 'FC2', 'C2', 'C4', 'C6', 'T8', 'CP4', 'P2', 'P8', 'P10', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'TP7', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F8', 'FT8', 'FC6', 'C4', 'C6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'TP7', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC2', 'TP8', 'CP6', 'P2', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'TP7', 'P1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F8', 'FCz', 'Cz', 'P2', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'T7', 'TP7', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'AF8', 'FC4', 'Cz', 'C2', 'C4', 'T8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['AF7', 'TP7', 'AF8']
    Rejecting  epoch based on EEG : ['P9', 'O1', 'Oz', 'AF8']
    Rejecting  epoch based on EEG : ['AF4', 'P10']
    Rejecting  epoch based on EEG : ['P9']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'O1', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'P10']
    Rejecting  epoch based on EEG : ['AF7', 'T7', 'TP7', 'P9', 'Iz', 'POz', 'CPz', 'AF8', 'F2', 'F4', 'FC6', 'FC4', 'C2', 'C4', 'T8', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P9', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'FC4', 'C4', 'P4', 'P10', 'O2']
    Rejecting  epoch based on EEG : ['P7', 'Fp2', 'P10']
    Rejecting  epoch based on EEG : ['PO7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F2', 'F4', 'FC4', 'FCz', 'C2', 'C4']
    Rejecting  epoch based on EEG : ['PO7', 'PO3', 'O1', 'Pz', 'P2', 'P4', 'P6', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['TP7', 'P9']
    Rejecting  epoch based on EEG : ['AF4']
    Rejecting  epoch based on EEG : ['P10']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'O2']
    Rejecting  epoch based on EEG : ['P9', 'O1']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P9', 'O1', 'Fpz', 'Fp2', 'AF8', 'AF4', 'P10']
    Rejecting  epoch based on EEG : ['P7']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P9', 'Iz', 'Fpz', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['P1', 'P9', 'Iz', 'O2']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'O1', 'Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'AF8', 'P6', 'P8', 'PO8', 'PO4']
    Rejecting  epoch based on EEG : ['P7']
    Rejecting  epoch based on EEG : ['AF3']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'P4', 'P6', 'P8', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Iz', 'Oz', 'CPz', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'Iz', 'Pz']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'P7', 'Iz', 'Oz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'FC2', 'FCz', 'P8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'P9', 'O1', 'Iz', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FC4', 'FC2', 'FCz', 'P8', 'P10', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['F1', 'P9', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F2']
    Rejecting  epoch based on EEG : ['P5', 'P7', 'P9', 'PO7', 'PO3', 'Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['P8']
    Rejecting  epoch based on EEG : ['TP7']
    Rejecting  epoch based on EEG : ['P1', 'Pz', 'AF8']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'PO3']
    Rejecting  epoch based on EEG : ['Fp1', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'POz', 'Pz', 'CPz', 'P2', 'PO4', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'POz', 'Pz', 'CPz', 'P2', 'Temp']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'P5', 'P7', 'P9', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FCz', 'P10']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P9', 'O1', 'Oz', 'POz', 'Pz', 'CP2', 'P2']
    Rejecting  epoch based on EEG : ['TP7', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'Fz']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'O1', 'Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['P9', 'PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'O1', 'P10']
    Rejecting  epoch based on EEG : ['P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F4', 'F6', 'F8', 'FC4', 'FC2', 'P10', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['P9', 'AF8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'FC3', 'FC1', 'P7', 'P9', 'PO3', 'O1', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4']
    Rejecting  epoch based on EEG : ['F1', 'FC1', 'Fp2', 'FC4', 'FCz']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'O1', 'POz', 'Fp2', 'FC4']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'AF4', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['PO3', 'P4', 'P6', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P9', 'PO7', 'O1', 'Iz', 'Oz', 'Fp2', 'AF8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['P2', 'PO4']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'PO3', 'Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['PO7', 'PO3', 'Oz', 'POz', 'Fpz', 'Fp2', 'AF4', 'AFz', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz', 'PO8']
    Rejecting  epoch based on EEG : ['P1', 'P7', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'Fp2', 'AF8', 'AF4', 'P8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'P3', 'PO7', 'PO3', 'O1', 'Temp']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['PO4']
    Rejecting  epoch based on EEG : ['PO7', 'C4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8', 'AF4', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8', 'F8']
    Rejecting  epoch based on EEG : ['P2', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8', 'PO4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8', 'P8', 'P10', 'PO8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fz']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp1', 'Temp']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'F8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['PO7', 'O1']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp1', 'Temp']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8']
    Rejecting  epoch based on EEG : ['F8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['FC2']
    Rejecting  epoch based on EEG : ['F7']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'P8', 'P10', 'PO8']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'P8', 'P10', 'PO8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8', 'F8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'F8']
    Rejecting  epoch based on EEG : ['FC2']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['F8']
    Rejecting  epoch based on EEG : ['Fpz', 'AF8']
    Rejecting  epoch based on EEG : ['TP7', 'P9', 'PO7', 'AF8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['POz', 'Fp2', 'AF8', 'P4', 'PO4']
    Rejecting  epoch based on EEG : ['Fp1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'Iz', 'Oz', 'Fpz', 'AF8', 'AF4', 'AFz', 'F2', 'F8', 'FT8', 'FC6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'AF4', 'FC4', 'CP2', 'P2', 'P4', 'P6', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'P1', 'P7', 'O1', 'Oz', 'POz', 'Pz', 'P2', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['O2']
    Rejecting  epoch based on EEG : ['P6', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['TP7', 'CP1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Fp2', 'AF8', 'AFz', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'P8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'TP7', 'P1', 'P3', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'P1', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['FCz']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['CP4', 'P6']
    Rejecting  epoch based on EEG : ['P8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F3', 'F7', 'AF8', 'F8', 'FT8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['P4', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P5', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F8']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['PO8', 'PO4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F2', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['O1', 'Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['F7', 'FT7', 'F8', 'FC6']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF3', 'Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['AF8', 'F8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F8']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F2', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['TP8']
    Rejecting  epoch based on EEG : ['P9', 'O1']
    Rejecting  epoch based on EEG : ['AF7', 'P9', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['P8']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P9']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'FC2']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['PO7', 'Iz', 'Oz']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'O1', 'Temp']
    Rejecting  epoch based on EEG : ['P2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['CP1']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['PO8', 'PO4']
    Rejecting  epoch based on EEG : ['F8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['P4']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz']
    Rejecting  epoch based on EEG : ['P2']
    Rejecting  epoch based on EEG : ['FT7', 'P10']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'F8']
    Rejecting  epoch based on EEG : ['P4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'FC2', 'FCz']
    Rejecting  epoch based on EEG : ['FC2', 'FCz']
    Rejecting  epoch based on EEG : ['PO4']
    Rejecting  epoch based on EEG : ['O2']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'FT7', 'P1', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F6', 'FC6', 'FC4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['FT7']
    Rejecting  epoch based on EEG : ['FC2']
252 bad epochs dropped


%matplotlib inline
epochs.plot();


%matplotlib inline
epochs['standard'].plot_image()

Not setting metadata
375 matching events found
No baseline correction applied
0 projection items activated
combining channels using "gfp"

[<Figure size 432x288 with 3 Axes>]


%matplotlib inline
epochs['oddball_with_reponse'].plot_image()

Not setting metadata
84 matching events found
No baseline correction applied
0 projection items activated
combining channels using "gfp"

[<Figure size 432x288 with 3 Axes>]


evoked_standard = epochs['standard'].average()
evoked_oddball_resp = epochs['oddball_with_reponse'].average()

NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).
NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).


mne.viz.plot_evoked(evoked_standard, spatial_colors = True, ylim = dict(eeg = [-20, 20]));


mne.viz.plot_evoked(evoked_oddball_resp, spatial_colors = True, ylim = dict(eeg = [-20, 20]));


evoked_standard.plot_joint();

No projector specified for this dataset. Please consider the method self.add_proj.
NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).


evoked_oddball_resp.plot_joint();

No projector specified for this dataset. Please consider the method self.add_proj.
NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).


mne.viz.plot_compare_evokeds([evoked_standard, evoked_oddball_resp], axes = 'topo');


def eeg_prepross(subj, run):
    # get subj number (e.g., 001, 002, 003...)
    sub_num = subj
    
    # Read the set file
    raw_path = f"./oddball/sub-{sub_num}/eeg"
    raw = mne.io.read_raw_eeglab(raw_path + f"/sub-{sub_num}_task-P300_run-{run}_eeg.set", preload = True)
    
    # Select only a subset of the channels (only EEG channels)
    raw_eeg = raw.copy()
    raw_eeg.pick_types(meg = False, eeg = True, eog = True, exclude = [])
    
    # Referencing using the T7, T8 channels as the referenced electrode
    raw_ref = raw_eeg.copy()
    raw_ref.set_eeg_reference(ref_channels = ['T7', 'T8'])
    
    # Using a band pass filter from 0.16-30 Hz
    raw_filtered = raw_ref.copy()
    raw_filtered = raw_filtered.filter(l_freq = 0.16, h_freq = 30)
    
    # Get events and event ids from annotations
    events, event_id = mne.events_from_annotations(raw_filtered)
    
    # Create epoch for the events
    # Epochs will start 200 ms before the stimulus onset and end 
    # 600 ms after stimulus onset
    tmin = -0.2
    tmax = 0.6
    baseline = (None, 0)
    
    # Rejecting epochs with peak-to-peak amplitude +/-100 µV
    reject_criteria = dict(
        eeg = 100e-6,  # 100 µV
    )  

    epochs = mne.Epochs(raw_filtered,
                        events = events,
                        event_id = event_id,
                        reject = reject_criteria,
                        tmin = tmin,
                        tmax = tmax,
                        baseline = baseline,
                        preload = True)
    
    evoked_standard = epochs['standard'].average()
    evoked_oddball_resp = epochs['oddball_with_reponse'].average() 
    
    return(evoked_standard, evoked_oddball_resp)


# Loop over all 13 subjects to do the pre-processing for all three runs

evoked_dict = {}
# {
#     '001': [run1, run2, run3],
#     '002': [run1, run2, run3], 
#     ...
#     '013': [run1, run2, run3]

# }

# run1[0] = standard
# run1[1] = oddball

sub_index = ['001', '002', '003', '004', '005', '006', '007', '008', '009', '010', '011', '012', '013']
runs = range(1, 4)

for subj in sub_index:
    subj_runs = []
    for run in runs:
        subj_runs.append(eeg_prepross(subj, run))
    evoked_dict[subj] = subj_runs

NOTE: pick_types() is a legacy function. New code should use inst.pick(...).
EEG channel type selected for re-referencing
Applying a custom ('EEG',) reference.
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 0.16 - 30 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 0.16
- Lower transition bandwidth: 0.16 Hz (-6 dB cutoff frequency: 0.08 Hz)
- Upper passband edge: 30.00 Hz
- Upper transition bandwidth: 7.50 Hz (-6 dB cutoff frequency: 33.75 Hz)
- Filter length: 5281 samples (20.629 s)

[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done   1 out of   1 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   2 out of   2 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   3 out of   3 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   4 out of   4 | elapsed:    0.0s remaining:    0.0s

Used Annotations descriptions: ['ignore', 'noise', 'noise_with_reponse', 'oddball', 'oddball_with_reponse', 'response', 'standard']
Not setting metadata
863 matching events found
Setting baseline interval to [-0.19921875, 0.0] s
Applying baseline correction (mode: mean)
0 projection items activated
Using data from preloaded Raw for 863 events and 206 original time points ...
    Rejecting  epoch based on EEG : ['P10']
    Rejecting  epoch based on EEG : ['P10']
    Rejecting  epoch based on EEG : ['FT7', 'P10']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FT7', 'FC5', 'FC3', 'FC1', 'C1', 'C3', 'C5', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC3', 'FC1', 'C1', 'C3', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC5', 'FC3', 'FC1', 'C1', 'C3', 'T7', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'T8', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC3', 'FC1', 'C1', 'C3', 'T7', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'T8', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'F3', 'F5', 'F7', 'FC1', 'C1', 'T7', 'TP7', 'CP5', 'CP3', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4', 'C6', 'T8', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'F5', 'T7', 'TP7', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC4', 'FC2', 'FCz', 'C2', 'C4', 'C6', 'T8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'T7', 'P9', 'O1', 'Iz', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F8', 'FT8', 'FC6', 'FC2', 'C2', 'C4', 'C6', 'T8', 'CP4', 'P2', 'P8', 'P10', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'TP7', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F8', 'FT8', 'FC6', 'C4', 'C6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'TP7', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FC2', 'TP8', 'CP6', 'P2', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'AF7', 'AF3', 'F1', 'TP7', 'P1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F8', 'FCz', 'Cz', 'P2', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'T7', 'TP7', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'AF8', 'FC4', 'Cz', 'C2', 'C4', 'T8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['AF7', 'TP7', 'AF8']
    Rejecting  epoch based on EEG : ['P9', 'O1', 'Oz', 'AF8']
    Rejecting  epoch based on EEG : ['AF4', 'P10']
    Rejecting  epoch based on EEG : ['P9']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'O1', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'P10']
    Rejecting  epoch based on EEG : ['AF7', 'T7', 'TP7', 'P9', 'Iz', 'POz', 'CPz', 'AF8', 'F2', 'F4', 'FC6', 'FC4', 'C2', 'C4', 'T8', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P9', 'Pz', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'FC4', 'C4', 'P4', 'P10', 'O2']
    Rejecting  epoch based on EEG : ['P7', 'Fp2', 'P10']
    Rejecting  epoch based on EEG : ['PO7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F2', 'F4', 'FC4', 'FCz', 'C2', 'C4']
    Rejecting  epoch based on EEG : ['PO7', 'PO3', 'O1', 'Pz', 'P2', 'P4', 'P6', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['TP7', 'P9']
    Rejecting  epoch based on EEG : ['AF4']
    Rejecting  epoch based on EEG : ['P10']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'O2']
    Rejecting  epoch based on EEG : ['P9', 'O1']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P9', 'O1', 'Fpz', 'Fp2', 'AF8', 'AF4', 'P10']
    Rejecting  epoch based on EEG : ['P7']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P9', 'Iz', 'Fpz', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['P1', 'P9', 'Iz', 'O2']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'O1', 'Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'AF8', 'P6', 'P8', 'PO8', 'PO4']
    Rejecting  epoch based on EEG : ['P7']
    Rejecting  epoch based on EEG : ['AF3']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'P4', 'P6', 'P8', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Iz', 'Oz', 'CPz', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'Iz', 'Pz']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'P7', 'Iz', 'Oz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'FC2', 'FCz', 'P8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'P9', 'O1', 'Iz', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FC4', 'FC2', 'FCz', 'P8', 'P10', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['F1', 'P9', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F2']
    Rejecting  epoch based on EEG : ['P5', 'P7', 'P9', 'PO7', 'PO3', 'Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['P8']
    Rejecting  epoch based on EEG : ['TP7']
    Rejecting  epoch based on EEG : ['P1', 'Pz', 'AF8']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'PO3']
    Rejecting  epoch based on EEG : ['Fp1', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'POz', 'Pz', 'CPz', 'P2', 'PO4', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'CP1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'POz', 'Pz', 'CPz', 'P2', 'Temp']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'P5', 'P7', 'P9', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'FCz', 'P10']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P9', 'O1', 'Oz', 'POz', 'Pz', 'CP2', 'P2']
    Rejecting  epoch based on EEG : ['TP7', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'Fz']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'O1', 'Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['P9', 'PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'O1', 'P10']
    Rejecting  epoch based on EEG : ['P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F4', 'F6', 'F8', 'FC4', 'FC2', 'P10', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['P9', 'AF8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'FC3', 'FC1', 'P7', 'P9', 'PO3', 'O1', 'CPz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'F6', 'FC4', 'FC2', 'FCz', 'Cz', 'C2', 'C4']
    Rejecting  epoch based on EEG : ['F1', 'FC1', 'Fp2', 'FC4', 'FCz']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'O1', 'POz', 'Fp2', 'FC4']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'AF4', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['PO3', 'P4', 'P6', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P9', 'PO7', 'O1', 'Iz', 'Oz', 'Fp2', 'AF8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['P2', 'PO4']
    Rejecting  epoch based on EEG : ['P7', 'PO7', 'PO3', 'Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['PO7', 'PO3', 'Oz', 'POz', 'Fpz', 'Fp2', 'AF4', 'AFz', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz', 'PO8']
    Rejecting  epoch based on EEG : ['P1', 'P7', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'Fp2', 'AF8', 'AF4', 'P8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'P3', 'PO7', 'PO3', 'O1', 'Temp']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['PO4']
    Rejecting  epoch based on EEG : ['PO7', 'C4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8', 'AF4', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8', 'F8']
    Rejecting  epoch based on EEG : ['P2', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8', 'PO4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8', 'P8', 'P10', 'PO8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fz']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp1', 'Temp']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'F8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['PO7', 'O1']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp1', 'Temp']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8']
    Rejecting  epoch based on EEG : ['F8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['FC2']
    Rejecting  epoch based on EEG : ['F7']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'P8', 'P10', 'PO8']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'P8', 'P10', 'PO8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8', 'F8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']

[Parallel(n_jobs=1)]: Done  65 out of  65 | elapsed:    0.4s finished

    Rejecting  epoch based on EEG : ['Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'F8']
    Rejecting  epoch based on EEG : ['FC2']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['F8']
    Rejecting  epoch based on EEG : ['Fpz', 'AF8']
    Rejecting  epoch based on EEG : ['TP7', 'P9', 'PO7', 'AF8', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['POz', 'Fp2', 'AF8', 'P4', 'PO4']
    Rejecting  epoch based on EEG : ['Fp1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'Iz', 'Oz', 'Fpz', 'AF8', 'AF4', 'AFz', 'F2', 'F8', 'FT8', 'FC6', 'TP8', 'CP6', 'CP4', 'CP2', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'P1', 'P3', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'CPz', 'AF4', 'FC4', 'CP2', 'P2', 'P4', 'P6', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'P1', 'P7', 'O1', 'Oz', 'POz', 'Pz', 'P2', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['O2']
    Rejecting  epoch based on EEG : ['P6', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['TP7', 'CP1', 'P5', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Fp2', 'AF8', 'AFz', 'F4', 'F6', 'F8', 'FT8', 'FC6', 'P8', 'P10', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'TP7', 'P1', 'P3', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'P2', 'P4', 'P6', 'P8', 'P10', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'P1', 'P7', 'P9', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'Pz', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8', 'PO8', 'PO4', 'O2', 'Temp']
    Rejecting  epoch based on EEG : ['FCz']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['CP4', 'P6']
    Rejecting  epoch based on EEG : ['P8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F3', 'F7', 'AF8', 'F8', 'FT8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['P4', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['P5', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F8']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['PO8', 'PO4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F2', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['O1', 'Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['F7', 'FT7', 'F8', 'FC6']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF3', 'Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['AF8', 'F8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'F8']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F2', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['TP8']
    Rejecting  epoch based on EEG : ['P9', 'O1']
    Rejecting  epoch based on EEG : ['AF7', 'P9', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['P8']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P9']
    Rejecting  epoch based on EEG : ['P7', 'P9', 'PO7', 'FC2']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['PO7', 'Iz', 'Oz']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F8', 'O2']
    Rejecting  epoch based on EEG : ['Fp1', 'Temp']
    Rejecting  epoch based on EEG : ['Fp1', 'O1', 'Temp']
    Rejecting  epoch based on EEG : ['P2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['CP1']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['PO8', 'PO4']
    Rejecting  epoch based on EEG : ['F8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['P4']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Pz']
    Rejecting  epoch based on EEG : ['P2']
    Rejecting  epoch based on EEG : ['FT7', 'P10']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'F8']
    Rejecting  epoch based on EEG : ['P4']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F1', 'F3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'Fz', 'F2', 'F4', 'FC2', 'FCz']
    Rejecting  epoch based on EEG : ['FC2', 'FCz']
    Rejecting  epoch based on EEG : ['PO4']
    Rejecting  epoch based on EEG : ['O2']
    Rejecting  epoch based on EEG : ['AF7']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'FT7', 'P1', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F6', 'FC6', 'FC4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['FT7']
    Rejecting  epoch based on EEG : ['FC2']
252 bad epochs dropped
NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).
NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).
NOTE: pick_types() is a legacy function. New code should use inst.pick(...).
EEG channel type selected for re-referencing
Applying a custom ('EEG',) reference.
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 0.16 - 30 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 0.16
- Lower transition bandwidth: 0.16 Hz (-6 dB cutoff frequency: 0.08 Hz)
- Upper passband edge: 30.00 Hz
- Upper transition bandwidth: 7.50 Hz (-6 dB cutoff frequency: 33.75 Hz)
- Filter length: 5281 samples (20.629 s)

[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done   1 out of   1 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   2 out of   2 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   3 out of   3 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   4 out of   4 | elapsed:    0.0s remaining:    0.0s

Used Annotations descriptions: ['noise', 'noise_with_reponse', 'oddball_with_reponse', 'response', 'standard']
Not setting metadata
862 matching events found
Setting baseline interval to [-0.19921875, 0.0] s
Applying baseline correction (mode: mean)
0 projection items activated
Using data from preloaded Raw for 862 events and 206 original time points ...
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['PO8', 'O2']
    Rejecting  epoch based on EEG : ['C2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['FC5']
    Rejecting  epoch based on EEG : ['FC5']
    Rejecting  epoch based on EEG : ['FC5']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['P1']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['PO7']
    Rejecting  epoch based on EEG : ['PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['PO8', 'O2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'AFz', 'F4', 'F6', 'F8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp1', 'PO7', 'PO3', 'O1', 'Iz', 'Oz', 'POz', 'P8', 'PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['FT7', 'P10']
    Rejecting  epoch based on EEG : ['FT7']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF3']
    Rejecting  epoch based on EEG : ['AF3']
    Rejecting  epoch based on EEG : ['Fp2']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['P2']
    Rejecting  epoch based on EEG : ['Pz']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'Fpz', 'Fp2', 'AF8', 'AF4', 'AFz']
    Rejecting  epoch based on EEG : ['POz']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['Fpz', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['Fp2']

[Parallel(n_jobs=1)]: Done  65 out of  65 | elapsed:    0.2s finished

    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['FC6']
    Rejecting  epoch based on EEG : ['AF4']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8', 'AF4', 'F8']
    Rejecting  epoch based on EEG : ['C2']
    Rejecting  epoch based on EEG : ['AF3']
    Rejecting  epoch based on EEG : ['PO8']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['C3']
    Rejecting  epoch based on EEG : ['F8']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2']
    Rejecting  epoch based on EEG : ['AF8', 'Cz']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['PO8', 'O2']
    Rejecting  epoch based on EEG : ['AF8', 'F8', 'FC2']
    Rejecting  epoch based on EEG : ['PO7', 'O1', 'PO8', 'O2']
    Rejecting  epoch based on EEG : ['Pz']
    Rejecting  epoch based on EEG : ['C3', 'Pz']
    Rejecting  epoch based on EEG : ['C3']
    Rejecting  epoch based on EEG : ['PO8', 'PO4', 'O2']
    Rejecting  epoch based on EEG : ['PO8', 'O2']
    Rejecting  epoch based on EEG : ['AF7', 'AF8']
    Rejecting  epoch based on EEG : ['AF8']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz']
    Rejecting  epoch based on EEG : ['Fpz', 'Fp2', 'AF8']
    Rejecting  epoch based on EEG : ['FC2']
    Rejecting  epoch based on EEG : ['AF7', 'Fpz', 'Fp2', 'AF8', 'AF4', 'FC2']
    Rejecting  epoch based on EEG : ['AF7', 'AF3', 'F5', 'F7', 'AF8', 'F8', 'FT8']
    Rejecting  epoch based on EEG : ['O1']
    Rejecting  epoch based on EEG : ['O1', 'Oz']
    Rejecting  epoch based on EEG : ['AF7', 'AF8', 'F8']
    Rejecting  epoch based on EEG : ['F7', 'FT7', 'FC5']
    Rejecting  epoch based on EEG : ['FT7']
    Rejecting  epoch based on EEG : ['FT7']
82 bad epochs dropped
NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).
NOTE: pick_channels() is a legacy function. New code should use inst.pick(...).
NOTE: pick_types() is a legacy function. New code should use inst.pick(...).
EEG channel type selected for re-referencing
Applying a custom ('EEG',) reference.
Filtering raw data in 1 contiguous segment
Setting up band-pass filter from 0.16 - 30 Hz

FIR filter parameters
---------------------
Designing a one-pass, zero-phase, non-causal bandpass filter:
- Windowed time-domain design (firwin) method
- Hamming window with 0.0194 passband ripple and 53 dB stopband attenuation
- Lower passband edge: 0.16
- Lower transition bandwidth: 0.16 Hz (-6 dB cutoff frequency: 0.08 Hz)
- Upper passband edge: 30.00 Hz
- Upper transition bandwidth: 7.50 Hz (-6 dB cutoff frequency: 33.75 Hz)
- Filter length: 5281 samples (20.629 s)

[Parallel(n_jobs=1)]: Using backend SequentialBackend with 1 concurrent workers.
[Parallel(n_jobs=1)]: Done   1 out of   1 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   2 out of   2 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   3 out of   3 | elapsed:    0.0s remaining:    0.0s
[Parallel(n_jobs=1)]: Done   4 out of   4 | elapsed:    0.0s remaining:    0.0s


# Loop over all subjects to combine the evoked objects across 3 runs

combined_evoked_dict = {}
# {
#     '001': [standard, oddball],
#     '002': [standard, oddball], 
#     ...
#     '013': [standard, oddball]

# }

for subj in sub_index:
    combined = []
    for run in runs:
        # combined standard (run 1-3)
        combined.append(mne.combine_evoked([evoked_dict[subj][0][0], evoked_dict[subj][1][0], evoked_dict[subj][2][0]], 'equal'))
        # combined oddball (run 1-3)
        combined.append(mne.combine_evoked([evoked_dict[subj][0][1], evoked_dict[subj][1][1], evoked_dict[subj][2][1]], 'equal'))
    combined_evoked_dict[subj] = combined


# Subject 8 standard condition
mne.viz.plot_evoked(combined_evoked_dict['008'][0], spatial_colors = True, ylim = dict(eeg = [-20, 20]));


# Subject 8 oddball condition
mne.viz.plot_evoked(combined_evoked_dict['008'][1], spatial_colors = True, ylim = dict(eeg = [-20, 20]));


# Exclude subject 2 due to missing values
sub_exclude002 = [x for x in sub_index if x != '002']

def plot_all_evoked():
    for subj in sub_exclude002:
        standard = combined_evoked_dict[subj][0]
        oddball = combined_evoked_dict[subj][1]
        
        print(f"Subject {subj} (Condition: standard)")
        mne.viz.plot_evoked(standard, spatial_colors = True, ylim = dict(eeg = [-20, 20]));
        
        print(f"Subject {subj} (Condition: oddball)")
        mne.viz.plot_evoked(oddball, spatial_colors = True, ylim = dict(eeg = [-20, 20]));
        
        print("--------------------------------------------------------------------\n\n")


plot_all_evoked()

Subject 001 (Condition: standard)

Subject 001 (Condition: oddball)

--------------------------------------------------------------------


Subject 003 (Condition: standard)

Subject 003 (Condition: oddball)

--------------------------------------------------------------------


Subject 004 (Condition: standard)

Subject 004 (Condition: oddball)

--------------------------------------------------------------------


Subject 005 (Condition: standard)

Subject 005 (Condition: oddball)

--------------------------------------------------------------------


Subject 006 (Condition: standard)

Subject 006 (Condition: oddball)

Measurement date	Unknown
Experimenter	Unknown
Participant	Unknown
Digitized points	68 points
Good channels	65 EEG, 11 misc, 2 Galvanic skin response, 1 RESP
Bad channels	None
EOG channels	Not available
ECG channels	Not available
Sampling frequency	256.00 Hz
Highpass	0.00 Hz
Lowpass	128.00 Hz
Filenames	sub-001_task-P300_run-1_eeg.set
Duration	00:12:38 (HH:MM:SS)

Measurement date	Unknown
Experimenter	Unknown
Participant	Unknown
Digitized points	68 points
Good channels	65 EEG
Bad channels	None
EOG channels	Not available
ECG channels	Not available
Sampling frequency	256.00 Hz
Highpass	0.00 Hz
Lowpass	128.00 Hz
Filenames	sub-001_task-P300_run-1_eeg.set
Duration	00:12:38 (HH:MM:SS)

Measurement date	Unknown
Experimenter	Unknown
Participant	Unknown
Digitized points	68 points
Good channels	65 EEG
Bad channels	None
EOG channels	Not available
ECG channels	Not available
Sampling frequency	256.00 Hz
Highpass	0.00 Hz
Lowpass	128.00 Hz
Filenames	sub-001_task-P300_run-1_eeg.set
Duration	00:12:38 (HH:MM:SS)

Number of events	611
Events	ignore: 0 noise: 84 noise_with_reponse: 1 oddball: 0 oddball_with_reponse: 84 response: 67 standard: 375
Time range	-0.199 – 0.602 s
Baseline	-0.199 – 0.000 s

EEG Prepocessing in Python¶

BrainHack 2023 (Taiwan Hub)¶

Open dataset (OpenNeuro)¶

EEG data from an auditory oddball task¶

Read the raw data¶

Select only a subset of the channels¶

Visualize the sensor locations¶

Referencing¶

Filtering¶

Compare before and after filtering¶

Zoom in¶

Extract events from annotations¶

Creating epochs for the events¶

Creating evoked data¶

Preprocessing function¶